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ON TRAVELING WAVES IN BEAMS 


By Robert W. Leonard and Bernard Budiansky 


SUMMARY 


The basic equations of Timoshenko for the motion of vibrating non- 
uniform beams, which allow for effects of transverse shear deformation 
and rotary inertia, are presented in several forms, including one in 
which the equations &re written in the directions of the characteristics. 
The propagation of discontinuities in moment and shear, as governed by 
these equations, is discussed. 


Numerical traveling-wave solutions are obtained for some elementary 
problems of finite uniform beams for which the propagation velocities 
of bending and shear discontinuities are taken to be equal. These solu- 
tions are compared with modal solutions of Timoshenko's equations end, 
in some cases, with exact closed solutions. 


INTRODUCTION 


The theoretical analysis of transient stresses in aircraft wings 
and fuselages subjected to impact loadings has generally been performed 
by means of a mode-superposition method that uses the natural modes of 
vibration predicted by the elementary engineering theory of beam bending. 
(See, for example, refs. 1 to 3.) For very sharp impact loadings, how- 
ever, this approach is known to have certain shortcomings. For sharp 
impacts of short duration, many modes are often required for a satis- 
factory degree of convergence (see, for example, ref. 4); in addition, 
the use of elementary beam theory in the calculation of the higher modes 
of vibration is inaccurate because of the neglect of, among other factors, 
the effects of transverse shear deformation and rotary inertia which 
become increasingly important for higher and higher modes (ref. 5). 


A classically recognized alternative to the modal method of calcu- 
lating transient stresses in elastic bodies is the traveling-wave approach, 
which seeks to trace directly the propagation of stresses through the 
body (ref. 3). Although the traveling-wave concept has been successfully 
used to treat such simple problems as longitudinal and torsional impact 
of rods, only recently have serious attempts been made to study the 
transient bending response of beams by this approach. Flugge (ref. 6) 


- ~ - NACA TN 287} 


was apparently the first to point out that elementary beam theory could 
not serve as an adequate basis for the traveling-wave approach since the 
elementary theory predicts that disturbances propagate with infinite 
velocity; he showed, however, that a traveling-wave theory could be |. 
constructed by modifying the elementary theory, as Timoshenko (ref. 7) 
did, to include first-order effects of transverse shear deformation 

and rotary inertia. On the basis of this more accurate theory, Fligge 
found that discontinuities in moment and shear travel along the beam 
with finite, and generally distinct, velocities. A similar analysis 

was carried out independently by Robinson (ref. 8) who, exploiting the 
hyperbolic nature of Timoshenko's equations, proposed the use of approxi- 
mate methods of solution and gave some numerical results for a particular 
example. Pfeiffer (ref. 9) also suggested the possibility of step-by- 
step solutions by the method of characteristics. In reference 10, Uflyand 
attempted an analytical solution of Timoshenko's equations for the case 
of a simply supported beam subjected to a sudden application of load; 
however, as was shown by Dengler and Goland (ref. 11), Uflyand's work is 
marred by the fact that he applied boundary conditions that are incorrect 
for Timoshenko's theory. The only known example of an exact traveling- 
wave solution based on Timoshenko's theory was presented by Dengler and 
Goland for the case of an infinitely long beam meet ras to a concentrated 
impulse. 


Thus, although the use of Timoshenko's theory as a basis for the 
trensient-stress analysis of beams has been seriously considered, few 
problems have actually been solved. Additional basic studies of 
Timoshenko's equations and their solution, particularly for finite- 
length beams, constitute necessary prerequisites to the development of 
practical methods of dynemic-stress analysis based on the traveling-wave 
approach, In the present paper, several specific problems of transient 
loading of uniform beams of finite length are considered and their solu- 
tions by various procedures, all based on the Timoshenko theory, are 
presented. These procedures are (a) numerical step-by-step integration - 
the "method of characteristics," (b) mode superposition, and (c) exact 
closed-form solution. The examples are all for the special case of equal 
propagation velocities of shear and bending disturbances; only for this 
case have exact solutions been found in closed form. For the sake of 
completeness, the presentation of these solutions is preceded by an 
exposition of the basic equations of Timoshenko's theory, a derivation 
of the characteristic lines and characteristic forms of these equations, 
and a discussion of their implications concerning propagation of 
disturbances. 


SYMBOLS 


A cross-sectional area 


E Young's modulus of elasticity 
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G shear modulus of elasticity 


I cross-sectional moment of inertia 

E length of finite beam (arbitrary length for infinite beam) 
M internal bending moment (see fig. 1) 

M dimensionless internal bending moment, ML/Elp 


is] 


V vertical shear force on & cross section (see fig. 1) 
V dimensionless vertical shear force, WS) EIB 

EI 
Cy propagation velocity of bending discontinuities, = 
Co propagation velocity of shear discontinuities, FE 
p Operator used in the Laplace transformation 
q intensity of distributed external loading 
q dimensionless intensity of external loading, aL) Elp 
x cross-sectional radius of gyration 
t time 
v velocity of deflection, yr 
v dimensionless velocity of deflection, v/co 
x coordinate along beam 
y deflection (see fig. 1) 
y dimensionless deflection, y/L 

denotes & discontinuity in quantity immediately following 

is S x 
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E dimensionless coordinate along beam, x/L 

p density of beam material 

T dimensionless time, cjt/L 

Y rotation of cross section about neutral axis (see fig. 1) 
ay velocity of rotation of cross section, Vit 

a) dimensionless velocity of rotation, aL/c 7 

Subscripts: 

B contributing to resistance of beam to bending 

1 contributing to inertia 

S contributing to resistance of beam to shearing 


x,t,5,T indicate partial derivatives with respect to those quantities 


BASIC EQUATIONS 


Nonuniform Beams 


Timoshenko's equations.- The Timoshenko theory of beam bending 
(ref. 7) constitutes a modification of elementary beam theory that 
attempts to account for the effects of transverse shear deformation and 
rotary inertia; the basic assumption of elementary theory - that plane - 
sections remain plane - is retained.. The moment M, shear V, deflec- 
tion y, and cross-sectional rotation y of a nonuniform beam subject 
to a dynamic lateral loading q are governed, according to this theory, 
by the following four simultaneous partial-differential equations (see 


fis, Lys 


M +-EI,v, = 0 | (1a) 
V - AgG(yx - Y) = 0 (1b) 
Mx - V + platy, = 0 (1c) 
Vx 7 Pytt + 4 = O (14) 
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The first two equations constitute elastic laws relating the deformations 
to the internal loading. Equation (la) expresses the same relationship 
between moment and cross-sectional rotation as that given by elementary 
beam theory. Equation (1b) stipulates & linear relationship between the 
shear V and the shear angle yx -ýy at the neutral axis; Ag is the 
so-called "effective" shear-carrying area, different from the total 

area A4 since the true shear angle actually varies over the cross 
section. Equations (1c) and (1d) prescribe rotational and translational 
equilibrium, respectively, with the term pIi¥y+~ representing the con- 


tribution of rotary inertia. 


The moment and shear may be eliminated from equations (1) to yield 
two simultaneous partial-differential equations in y and Y: 


(Elpvrx), i AgG( Ya =) s els Yi = O 
(2) 


Aane - wm]. - PÁjYig +2 = 0 


This form is convenient for finding the normal modes and frequencies of 
free vibration (q = 0) predicted by Timoshenko's theory and for carrying 
out modal analyses that make use of these modes. In this theory, 4 


natural mode is described by a pair of functions [y(x),¥(x)] rather 
than a single function y(x), as is the case in elementary beam theory. 


For a traveling-wave analysis, however, it is advantageous to return 
to the original system of equations (eqs. (1)) but to put them in a more 
convenient form by differentiating equations (la) and (1b) with respect 
to time. The equations become: 


Oy + zz My = 0 (3a) 
Vy - E Vt -020 ( 3b) 
My. + pIj - Y = 0 (3c) 
V. - PAjvy +4 = 0 (3d) 


where the new variables, linear velocity v and angular velocity 4, 
have been introduced for Y. and Ves respectively. 
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Equations (3) comprise four first-order linear partial-differential 
equations in the four variables v, wœ, M, and V. Furthermore, equa- 
tions (3a) and (3c) contain derivatives of only M and 0, whereas 
equations (3b) and (3d) contain derivatives of only v and V. These 
facts are exploited in the next section in seeking characteristic lines, 
and characteristic forms, of these equations. 


Characteristics and the characteristic form of the-equations.- In 
equations (3a) and (3c), the variables M and o are differentiated 
with respect to both space and time; it would be advantageous to replace 
them by equivalent equations each involving only total derivatives (or 
differentials) in a particular direction in the space-time plane. The 
lines in the space-time plane having these particular directions - the 
characteristic lines or so-called "characteristics" - and the equivalent 
equations written in these directions are found as follows (ref. 12). 


A linear combination of equations (3a) and (3c), 


LR; + Elgox) + M. + ply - V=O (4) 


is formed, where the function u is to be determined in such a way 
that the partial derivatives in equation (4) combine to give total 
day 


derivatives = and do in the direction of an as yet unknown charac- 


teristic line [x(0),t(0) . In order that the terms involving deriva- 
tives of M combine in the form 


dx at dM 
Mr ac * "b gg " dg 


the function u must satisfy the following equation: 


where EM is the required slope of the characteristic line. Similarly, 
in order that the terms involving derivatives of œ combine in the form 


dx ¿a 3t _ dw 
OX do ' Ot dc — ao 
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the function uw must also satisfy the equation 


Since the characteristic slope must be the same in both cases, u is 
defined by 


ol; 
p° ~ Elp 
B 


Thus, the two values of yu and the corresponding characteristic 
slopes o are: 


dx 
H oe 
el 
(5) 
dt i1 
dX c4 
pa 
vl 
(6) 
du di 
dx Ci 


T - 

where c] = Si Then, multiplying equation (4) by dt and using equa- 
i 

tions (5) ylelds 


= aM + pI, do - V dt - 0 
"i 


when 
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Similarly, using equation (6) gives 


" = aM + oly do - V dt=0 
L 


when 


In an analogous fashion, it can be show, from equations (3b) 
and (3d), that - 


A 
1 aS . 
do dV - pAj dv + (a + E da = 0 
when 
at = - ax 
Co l 
and 
1 AgG 
- — dV - pA, dv + - — ojdt^ 0 
Co pr av Co ) 
when 
bs ea 
co 


AgG 
where Co = Era 
i 


Then, the system of equations (3) has associated with it four. real 
characteristic directions and is thus "totally hyperbolic" (ref. 12). 
A network of characteristics.can readily be constructed without prior 
knowledge of the unknowns M, V, 0, and v since their slopes depend 
only on the material and geometrical properties of the beam. For uniform 
beams, as well as for tapered beams having uniform material properties 
and geometrically similar cross sections, cj and co are constant, 


and the characteristics will therefore be straight; in general, however, 
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the characteristics may be curved. Figure 2 illustrates the four char- 
acteristics passing through a point P in the space-time plane with the 


characteristics designated as 


T at 1 
dx C1 

I-: at = = EJ 
dx Ci 
dx Co 
át l 

mons dx ^ ^ c5 


It is known that, by virtue of the totally hyperbolic character of the 


basic equations, the values of the unknowns M, V, 


wm, 


point P depend only on their initial values at t=0 between the 
points x; and xo on the beam (ref. 12). Furthermore, these values 


at P can, in turn, have influence only on points lying in the region 
above P enveloped by the I+ and I- characteristics through P. 
Thus no signal can proceed along the beam with a velocity greater than 
(which is generally larger than co, as illustrated in fig. 2). 


and v at the 


For the sake of easy reference, the four characteristic differential 


forms of the basic equations are grouped below. 


Along I+: + aM + pI, do - V at =0 
s1 
1 
Along I-: — aM - pl; do + V dt =0 
1 
Along II+: pe oA; dv + (pA, cao + a) dt = 
Al II-: 2 Gy oA pue VR dt = 


Related forms of these characteristic equations have 
Robinson (ref. 8) and Pfeiffer (ref. 9). 


been written by 


(Ta) 


(Tb) 


(7c) 


(Ta) 


=I 
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Propagation of discontinuities.- Characteristics are lines across 
which discontinuities may exist (ref. 12); indeed, this property 18 
often taken as the basic definition of a characteristic. In the present 
problem, discontinuities (or jumps) in M and «c can therefore exist 
across the I+ and I- characteristics, and discontinuities in V 
and v can exist across the II+ and II- characteristics. Hence, 

a jump in M or o will propagate with the velocity cy, whereas a 
Jump in V or v must proceed with the velocity co. It should be 
remarked that such discontinuities would appear only in the limiting 
cage of & beam subjected to an instantaneous loading. The solution of 
such idealized problems, which are often instructive, requires a knowl- 
edge of the laws governing the variations in strength of these discon- 
tinuities as they propagate through the beam. These laws are determined 
below for nonuniform beams for which it is assumed that the condition 

C4 = Cp does not hold over any portion of the beam; in other words, 


the I and II characteristics are distinct. The special case where 
Cj = Cp is considered subsequently when uniform beams are discussed. 


Let a and b be two points on a I- characteristic on either 
side of a particular I+ characteristic 


If M is discontinuous across the I+ characteristic, then M, - My 


retains a finite value OM as a and b are allowed to approach the 
I+ characteristic from either side. Consequently, from equation (7b) 
written along the I- characteristic, 


E 


since dt approaches zero as a and b approach each other. Thus, 
everywhere along a I+ characteristic, jumps 9M and 5w across this 
characteristic are related by 


SM = cioI, Sn | (8) 
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Similarly, the jumps across the eS characteristics can be readily 
shown to satisfy 


Along I-: BM = -cypIj dm (9) 
Along II+: SV = -copA,. Sv | (10) 
Along  II-: SV = c pA; dv (11) 


A jump in M is thus always accompanied by a definite jump in ao; 
similarly, jumps in V and v are always coupled together. 


The variations in the magnitude of & discontinuity as one proceeds 
along & characteristic may be determined in the following manner.  Equa- 
tion (7a) is written for the upper side and the lower side of the 
I+ characteristic; then, since V is continuous across a I+ char- 
acteristic, the difference of the two equations yields 


d(5M) T cy ply d( Sw) = 


along I+. Eliminating Sw by using equation (8) gives 


a(8M) + cor, al] = o 
cipi; 


By carrying out the indicated differentiation in the second term and 
dividing by 25M, the following result is obtained: 


a(sM) _ 1.4(c11;) 


S Meal So 


Solution of this equation gives 


(8M), = (8M)7 (12) 


as the relationship between the magnitude of the jumps in M at two 
points 1 and 2 onthe I+ characteristic. 
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It can be shown that the identical relationship holds between 
jumps in M at two points on a I- characteristic. Similarly, it can 
"be found that, on II+ and II- characteristics, 


(V), = (SV), (13) 


for any two points 1 and 2 on a given characteristic. The corre- 
sponding variations of the jumps Sw» and Ov are, of course, readily 
determined from equations (8) to (11). 


Uniform Beams 


Nondimensional form of the equations.- The examples to be presented 
in this report are all concerned with beams having uniform cross-sectional 


size and shape and uniform material properties throughout their length. 
For such beams, it is convenient to express Timoshenko's equations in 
nondimensional form. Thus, equations (3) may be written in terms of 


M = gr» 7- HE, B= TE, and ye as = 
Be +M, = 0 (15a) 
CY -RY-8-0 0007 -— 
Me +@,-V=0 (1hc) 
Vy -R = 7, +80 (1ha) 7 
where LE T "ed 1E. and TM The quantity L 


refers to the beam length for all beams except those of infinite length, 
in which case any convenient arbitrary length may be chosen for L. 


The characteristics of Timoshenko's equations for a uniform beam "E 
are defined in the E,r plane by the families of straight lines 
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I. EG =1 (15a) 
I-: i: sl (15b) 
TI+: it = = S (15c) 
II-: oe = | | (158) 


The nondimensionalized characteristic forms of.the basic equations become 


Along I+: dM + dx - VáT=0 , (168) 
Along I-: | dM - du +Y dr < 0 (16b) 
-— E Co 
Along II+: dv - R dv + (Fas + ES ajar = 0 (16c) 
ai = -- CA _ 
Along  I1-: dV + R dv + & Pm ajar = 0 (16d) 
E 


In addition to the. restriction to uniform beams, for which cy 
and co are constants, the examples presented herein are further limited 
to those beams for which the propagation velocities cj and cy are 
equal. This assumption has been made because the simplifications that 
result not only permit the ready attainment of numerical. solutions but 
also, in particular cases, permit the attainment of exact closed solu- 
tions for comparison. Since, for this very special case, the charac- 
teristics II coincide with the characteristics I, equations (16) may 
now be written, for q = O, as 


dp - Y at =0 
Along I+: : 3 (17a) 
av - LA a7 + IR aT 


S 


il 
O 


dM - db L YV dT = 0 
Along I-: (17b) 
av + hAPayv + LAD dr 


il 
O 
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Propagation of discontinuities when cj - co.- Equations (8) to (13), 

. which describe the behavior of discontinuities in a nonuniform beam, also 

describe, as a special case, the behavior of discontinuities in a uniform 
beam for which cy f cp. They show that such discontinuities propagate 


with constant magnitude. 


However, when the beam has properties such that cy = co, these 
equations are no longer valid. Equations which are valid in this case 
may be derived in precisely the same way by using equations (17) instead 
of equations (7). Discontinuities in such beams can be shown to be 
related by 


5M = 585 

Along I+: B (18a) 
BV = -hA?8y 
SM = -50 

Along I-: ` (18b) 
BV = A oF 


as they propagate, and they can be shown to vary in magnitude according 
to the equations 


duc Cb as 
LRR) - 5 V =0 
Along I+ and I-: . (19) 


d rar Dar _ 
a; (V) + 24780 = O 


Equations (19) may be solved simultaneously to obtain 


Su = A cos(XT - B) 
Along I+ and I-: (20) 


SV = -224 sin(XT - B) 


where A and B are arbitrary constants which must “be evaluated by 
using known values of 50 and SY at some point on the characteristic. 
The variations in ôM and SY can then be readily found by using 
equations (18). 


3B 
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Thus, for the case Cy = co, discontinuities in a uniform beam do 


not propagate unchanged but vary in magnitude sinusoidally as they 
progress through the beam. 


Limitations of the Theory 


It may be well to insert a word of caution about the applicability 
of Timoshenko's theory. The investigations of Prescott (ref. 13) and 
Cooper (ref. 14) have shown that, when the response of a beam includes 
components whose wave length is small compared to the depth of the beam, 
the assumption of Timoshenko's theory that plane sections remain plane 
after bending is, as might be expected, too restrictive. Since applied 
disturbances which could give rise to discontinuities would obviously 
excite even the shortest wave length in the spectrum of the response, 
the results obtained by application of Timoshenko's theory to such 
hypothetical problems cannot, in themselves, have practical significance. 


However, carrying out solutions involving discontinuities is a useful means 


of testing methods of solution of the Timoshenko equations with a view to 
applying these methods to problems in which discontinuities do not exist. 
Furthermore, the admittedly inaccurate response to an infinitely abrupt 
disturbance may be used to obtain the correct response to disturbances 

of a more practical nature through the application of Duhamel's super- 
position integral. 


SPECIFIC EXAMPLES - FINITE UNIFORM BEAMS WITH cy, = Co 


Three specific examples are considered; they are: a cantilever 
beam given a step velocity at the root, a simply supported beam subjected 
to a step moment at one end, and a simply supported beam subjected to a 
ramp-platform moment at one end. 


Methods of Solution 


In the examples, the results of calculation by the following three 
methods are compared: (a) numerical step-by-step integration along the 
characteristics, (b) normal-mode superposition, and (c) exact closed- 
form solution. The first two procedures are approximate in character, 
but they could conceivably be generalized sufficiently to be applied to 
practical structures; the last procedure, although exact, would rarely 
be useful in practice and is introduced herein chiefly as a check on the 
accuracy of the first two. 


The detailed descriptions and applications of the three methods are 
contained in appendixes A, B, and C. In brief, the numerical procedure 
exploits a grid of characteristic lines as shown in figure 3(a). For 


the case cj = Co that is under consideration, this grid consists of 


two families of lines in the £ r plane with slopes 1 and -l. A 
recurrence formula is developed in appendix A that gives the values 

of © and V at station 1 of a typical interior mesh (see fig. 3(b)) 
in terms of the values of © and V at stations 2, 3, and h. Repeated 
application of this formula, together with the use of special formulss 
for the half-meshes at either end of the beam and the knowledge of the 
magnitudes of jumps in O and V across characteristics where they 
occur, permits @ and V to be calculated throughout the ¢,7r plane. 
Subsequent determination of M and Y is achieved by means of simple 
addition formulas utilizing these calculated values of © and V. 


Although the solutions derived in appendix B have actually been 
obtained by Laplace transform techniques, they have been termed "modal 
solutions" because they are exactly those which would result from an 
application of the usual mode-superposition process. The exact closed 
solutions in appendix C have also been obtained through the use of 
Laplace transforms. 


Cantilever Beam Given a Step Velocity at the Root 


The first example to be considered, the response of a uniform 
cantilever beam given a step velocity v= 1 at the root, is the 
equivalent of the most severe idealized landing problem, the instan- 
taneous arrest of the root of a moving cantilever beam. Computed results 
for the shear and moment at the root-of such a beam having a ratio of 
length to radius of gyration of 10 (X = 5) and properties such that 
Cj = co are presented in figures b(a) and h(b), respectively. Two 


separate curves obtained by the numerical procedure are shown - one from 
a grid that divides the beam into 10 segments and the other from a 
20-segment solution. The modal solution includes the contributions of 
the first eight modes. Results given by an exact closed solution are 
shown for both the shear and moment at the root up to the time. T - 2. 
These exact results are actually those for an infinitely long beam, since 
the influence of the free end is not felt at the root until T= 2. 
After T= 2, the influence of the free end is felt and, in this case, 
an exact solution is not feasible. To illustrate the time range covered 
in the plots, the point corresponding to half the period of the first 
mode of vibration of the beam is indicated on the time scale of each 
plot. m 


In figure 4(a), the shear discontinuities evident in the numerical 
solutions occur each time the discontinuous wave front returns to the 
root after being reflected at the free end. The beam is seen to react 
violently to each of these boosts by the wave front, with more and more 
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oscillations occurring after each succeeding jump. The frequency of 
these oscillations tends to increase with each succeeding boost until 
limited by the finite time interval. In these regions of violent 
oscillation the &ccuracy of the numerical solutions is obviously ques- 
tionable; indeed the question arises &s to whether these oscillations 
are really predicted by the theory or are merely the result of some 
instability in the numerical process, This question is resolved in the 
next section in which the simply supported beam is considered. At any 
r&te, away from the regions of violent oscillations, the comparisons 
with the modal solution are favorable, and, for T < 2, the fine-grid 
numerical solution almost coincides with the exact closed solution valid 
in this region. 


The comparisons between the numerical and modal solutions are very 
good in figure (b) where the time history of the moment at the root is 
plotted. Again, the fine-grid numerical solution nearly coincides with 
the exact closed solution for T «2. 


Simply Supported Beam With an Applied End Moment 


Step moment.- A simply supported uniform beam with a ratio of 
length to radius of gyration of 10 (A = 5) and properties such that 
c] = co is subjected to a step moment M = 1 at the end $ - O0. Com- 


puted results for the shear at 5 - 0 are presented in figure 5(a) and 
the time history of moment at the center of the beam is presented in 
figure 5(b). The point corresponding to the full period of the first 
mode 1s marked on the time scale of each plot. The numerical curves 
for both shear and moment were obtained by a 20-segment solution. The 
modal curves were obtained by adding dynamic corrections to the static 
solutions, the dynamic corrections being calculated with six modes 

for both the shear and the moment. 


This example affords an answer to the question raised in the 
preceding section with regard to the stability of the numerical procedure 
following the passage of a discontinuous wave front. The fact is that 
the violent oscillations that occur after the discontinuity actually 
appear in the exact solution (fig. 5(a)) and are hence inherently present 
in the theory. | 


In figure 5(a) the scale of dimensionless vertical shear happens 
to be precisely the dynamic overshoot factor; it is of interest to note 
that values at least 15 times the static shear are predicted when the 
moment is applied suddenly. 


For shear (fig. 5(a)), the inaccuracies in the numerical solution 
Just after the discontinuities are evident; however, the numerical 
results approximate the exact solution very well elsewhere. A similar 
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Observation may be made for the modal solution, which, &s would be 
expected, ignores the discontinuities and violent oscillations caused 
by the wave front. The numerical and modal curves for the moment in 
figure 5(b) follow the same pattern, agreeing well everywhere except 
in the regions immediately following the discontinuities. 


Ramp-platform moment.- Computed results for the shear at the end 
t =0 and the moment at the center of the same simply supported uniform 
beam subjected to the applied ramp-platform moment 


M(O,T) = T | (OS 1 $ 1) 
M(0,T) = 1 (T > 1) 


are presented in figures 6(a) and 6(b). Again, the period of the first 
mode of vibration is marked on the time scales. 


The numerical curves were recomputed with a 20-segment grid for 
the new forcing function; however, the results for the modal and exact 
solutions were obtained by means of & superposition of the preceding 
results. This superposition was carried out analytically for the com- 
plete modal solution and for the exact solution in the range T <2. 
In the range T » 2, it was necessary to carry out the superposition 
for the exact solution numerically. B 


In figure 6 the time to peak value of the applied moment is 
seen to be approximately one-seventh the period of the first mode; 
predicted shear values approximately three and one-half times the 
static response occur. 


With the removal of the discontinuity, there remain no high- 
frequency oscillations which the numerical solution might be unable 
to represent. In fact, all three solutions for the shear and both 
moment solutions are seen to be in excellent agreement. 


CONCLUDING REMARKS 


Timoshenko's equations for the motion of vibrating nonuniform beams 
may be written in a characteristic form which appears to be well-suited 
to solution by numerical methods. In the examples presented in this 
report, all of which are for uniform beams with equal propagation veloc- 
ities of bending and shear discontinuities (cj = co), the solutions by the 
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numerical and modal methods generally agree well with each other as well 
as with exact closed solutions where these have been obtained. However, 
the modal method, of course, fails entirely to reproduce discontinuities 
in shear or moment, and the numerical procedure, although it yields these 
discontinuities, does so with decreasing accuracy as more and more 
reflections of the wave front occur. In the more practical situations 
where discontinuities do not exist, these difficulties will, of course, 


not arise. 


The results of the examples carried out by the numerical traveling- 
wave procedure encourage the viewpoint that such traveling-wave analyses 
may eventually be of practical usefulness. This kind of procedure is 
inherently simple and straightforward and has the advantage that the 
bulk of the labor involved is routine computation rather than math- 
ematical analysis. It should be emphasized, however, that numerical 
solutions of Timoshenko's equations have been demonstrated only for 
uniform beams in which the propagation velocities cj and co are 
equal; numerical procedures for the general case where they are unequal 
remain to be developed and tested. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., October 28, 1952. 
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APPENDIX A 


NUMERICAL SOLUTIONS FOR UNIFORM BEAMS WITH C] = Co 
Matrix Formulation 


Let the differential equations (17) be replaced by the finite- 
difference equations 


AM + AD - VAT = 0 (Ala) 
Along I+: = o 

AV - UÉAT + JAR Ar = 0 (Alb) 

MM - AB + V At = 0 (Alc) 
Along I-: 

AV + MEAT + hà S Ar = 0 (Ald) 


and consider a closely spaced network of I+ and I- characteristics - 
in the space-time plane as shown in figure 3(a). Let the corners of a 

typical interior mesh of this grid be designated as shown in figure 3(b). 

A step-by-step integration formula for œ and V may now be derived in = 
the following manner. D. = 


Equations (Al) may be written arene the upper sides of the typical 
mesh to obtain 


——— (A2) 
M - M -5 +0 + L(V + V3) =0 (A3) 
Vi - Vo - la^ (94 - Ve) + 22%AT(D] + DB) = 0 (A) 
V - V3 + la^ (v, - 73) + SPALE + 83) = 0 (A5) 


where © and Y have been assumed to vary linearly in the small 
intervals between the corners. Obviously equations (A2) to (45) | may _ . 
be solved simultaneously to obtain the four quantities V, O, V, and M 
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at point 1 in terms of their values at paints 2 and 3. However, it is 
noted that, if V, UD, V, and M at points 2 and 3 were determined by 
& similar process for the preceding meshes, they already satisfy the 
equations 


M, - My - Bp +0, + Vo + Vy) = 0 (A6) 
Mz - My + Bz - Dy - £z (v4 + vi) = 0 (AT) 
Vo - Vi + la (vo - Vi) + SPAR, + 3) = 0 (A8) 
Va - Vy - hA (v. - Vy) + 2A Ar (i. +y) = 0 (49) 
Equations (A6) and (AT) may be added to equations (A2) and (A3), 
respectively, to obtain 
a im i M = At i= = 
M -M +0, - +o, -F - V) = 0 
and 
— — E S MEN L d — 
M, - My, -© + 2D; - y + S(¥, - Vi) = 0 
which may, in turn, be subtracted to obtain the single equation 
B, - By - G, +B, - uv - WU 20 (A10) 


Similarly, equations (A8) and (A9) may be subtracted from equations (Ak) 
and (A5), respectively, to obtain 


ll 
O 


Vi - Wo + Vy TATE, - y + 2DMLAT(D, - 0 
h 1 à (à 4 


Vi - Wa + Vy + A (F - vy) + SAL - 23) 


I 
O 
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and these may be added to derive 
Z T7 v.c S O 
Vi - Vo - V3 + Vy + SANL - 3) = 0 (A11) 


Equations (ALO) and (All) may now be solved simultaneously to obtain 
the step-by-step integration formula 


5 : mod By 
NE Lr xe. M e ALI. (A12) 
Vz] (A Ar) +1 Va + V. Vi 
where 
AT 
1 2 
les} 7| 
“22, AT 1 
and 
L (X Ary? - 1 -AT 
Ay, | = 
LA AT (X Ar) - 1 


Thus ®© and V may be determined at every interior mesh point 
in the space-time plane by the repeated application of formula (A12) to 
each mesh as it is encountered. The half-meshes which occur at the 
vertical left and right boundaries of the plane (fig. 3(a)) require 
special formulas incorporating the known boundary data. These formulas 
may be derived from equations (Al) by a procedure similar to that used 
in obtaining equation (A12). Boundary formulas for some specific 
examples are presented in subsequent sections of this appendix. 


Besides boundary data, initial data must be provided in order that 
a step-by-step solution may be begun. In all the examples considered, 
the beam is initially at rest and has a disturbance applied at the 
point £=0 beginning at time t= 0. The region T < & (without grid 
lines in fig. 3(a)) therefore is one of zero stress and motion and Y 
and V are given along the line T = & by the known conditions at the 
wave front. It is with these values that each numerical solution is 


begun. 
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Discontinuities in ® and V offer no special difficulties since 
they propagate so as to be always located on I characteristics which, 
Of course, can be made part of the basic grid. Thus they are simply 
added as they are encountered. 


Once U and V have been determined at the grid intersections by 
the step-by-step process, M and Y may be found by repeated applica- 
tion of equation (A2) or (A3) and equation (Ah) or (A5) along the proper 
grid lines from boundaries where M and v are known. For example, 
if M is known along the left boundary, it may be found at a point g 


by applying equation (A2) successively to the intervals ab, bc, 
cd, . . . fg. The resulting expression for Mo becomes 


Mg = Ma + Gy - Gy + AS Va + Vp... + Ve +3 Vg (A13) 


This procedure is seen to correspond to integration of the first of 
differential equations (17a) by means of the trapezoidal rule. 


Specific Problems 


Cantilever beam given a step velocity at the root.- If the 


root E=0 of a uniform cantilever beam is given a step velocity v=1 
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at time T = O, the boundary conditions may be written 


v(0,T) = 1 
w(O0,T) = 0 
(A14) 
M(1,T) =0 
V(1,T) = 0 


Typical left- and right-boundary half-meshes are 


v(O, 7) 
WCO, Y) 


Application of equations (Alb) and (Ald) to the sides of the left- 
boundary half-mesh and proper combination of the resulting equations 
to eliminate "v produces 


Vi = 2V3 - Vh (A15) 


Similarly, application of-equations (Ala) and (Alc) to a right-boundary 
half-mesh produces 


so that a complete set of integration formulas for determining © and V 
for this problem is now available. In addition, if- the unused equations 
(eqs. (Ala) and (Alc) for the left-boundary half-mesh and eqs. (Alb) 

and (Ald) for the right) are applied and combined to eliminate M and Y 
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at points 3 and 2, the recurrence relations 

= = _ AT; s = 

M, = My - 25, - SV] - Vi) (A17) 
for left-boundary half-meshes and 


a xli. = = 
Y] = Y, - L8 Vo + Sa - Gy) (A18) 


for right-boundary half-meshes are obtained. These equations may be 
used to compute M at the root and V at the tip after © and V have 
been determined everywhere. 

Since the applied disturbance is eai d discontinuous, 0 and V 
will be discontinuous along the line = E (see fig. 3(a)). In order 
to begin the step-by-step solution fog. ( and V, these discontinuities 
must be determined in advance. Furthermore, since the wave front is 
reflected back into the beam at either end, © and V will also be 
discontinuous along the lines T = 2 - E, =2+ È... and O0 
and SV must be determined along each of uci lines before the step- 
by-step solution can be extended beyond it. The discontinuities at the 
wave front are determined as follows. 


From the boundary conditions (Alb), it is seen that 85(0,0) = O 


and 8v(0,0) = 1. Thus, from equation (18a), 68V(0,0) = -h3?, and 
equations (20) become, for the line T = E, 


&p(&£,t) = -2X sin AE 
e o (A19) 
SV(&,€) = -hA^cos AE 

At €=1, the discontinuities across = E are 6a(1,1) = 8M(1,1) = 


-2. sin X and 5V(1,1) = -4AF cos X, SO dixe if the boundary condi- 
tions M(1,T) = V(1, T) = 0 are to be satisfied, jumps must occur across 
=2- E at the point (1,1) with the magnitudes (1,1) = -8M(1,1) = 
-9X sin X and 5V(1,1) = hA^cos X. In view of these initial condi- 

tions at T = 1l, equations (20) become, for the line T - 2 - E, 


SRL £,2- E) 


5V(£,2-£) 


-2X sin AE 


(A20) 


hA*cos rE 
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Initial jump values . &9(0,2) = O and 8¥(0,2) = A? for the line 

T= 2+ E may be found from equation (A20) by satisfying the boundary 
conditions 0O(0,T) = 0 and v(0,T) =1. From these initial values, 
the discontinuities across T2 2-4 E are found to be the negative of 
those across T= &, or 


Sw(E,2+€) = 2X sin AE 
= (A21) 
BV(t,24&) = 4Afcos XE 
Then, it must be true that 
Sw( £, 4-6) = 2X sin AE 
(A22) 


Ecos Aj 


5v( E, h- E) 


and so forth, with the values on each succeeding line T =n * E repeating 
the values on the line T-(n- lD TL t, The variations in the magnitudes 
of BD and BV as the wave front propagates back and forth through the 
beam are thus clearly defined, and, since T <L is a region of zero 


stress and motion, equations (A19) define the values of œw and V along 
the line T= £, 


With © and V known along T= E, the solution for d and V 
may be begun by applying formula (Al5) to the half-mesh in the lower 
corner of triangular region (1) (fig. 3(a)). It is continued step by 
step throughout the triangle, with formula (A12) being used for all 
interior meshes. When the first triangle is complete, the known jumps 
along the line vr = 2 ~- E may be added to the values computed for the. 
under side of this line, and the solution may then be carried out in 


triengle (2) beginning again at the lower corner this time with 
formula (A16). In this way the solution may be carried through as many 
triangular regions as desired. 


Two sets of computations, with time intervals AT = 0.1 and 
AT = 0.05, have been made for a cantilever beam for which X= 5. These 
time intervals correspond to grids dividing the beam into- 10 and 
20 segments, respectively. (The computations for the 20-segment solution 
were made on the Bell Telephone Laboratories X-66744 relay computer at 
the Langley Laboratory.) The computed time histories of shear and moment 
at the root have been plotted in figures h(a) and 4(b), respectively, up 
to t= lO. In order to obtain these time histories, the computations 
for @® and V had to be carried through the first nine triangular 
regions. l l 
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Simply supported beam with an applied end moment.- If the end E = 0 
of a uniform simply supported beam for which cj = co is subjected to 
an applied bending moment Mo(T) beginning at time T = 0, the boundary 
conditions are 


M(0,T) = Mo (tT ) 

v(O,T) = O 

- (A23) 
M(1,T) = 0 

v(1,r) =0 


Typical left- and right-boundary half-meshes for this beam are 


M(O, T) 
v(0, 7) 


fa = 
sn 9,97 
G H 
Se e 
d 
TN 
o O 


Equations (Alc), (Ala), (Ald), and (Alb) may be written &long the 
sides of the left half-mesh to give, respectively, 


(Mo), - Ma - 3 +5 + NV, + V3) = o 
Mz - (Mo), +53 - 5, - S{V3 + Ty) = 0 


I 
O 


| 
O 


¥3 - Vj - Pv + SED, + Dy) = 
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Adding the first two equations and ERTAN 1 the last two equations 
produces um 


(Mo); - (Mo), - B + 253 - à, + FV - vi) = o 


and 
V, - EV. + Vy + SAE, - 3) = 0 


which may be solved simultaneously for < and Vj to obtain 


0 1 20) A 0) " = l h 
Va, j (A Ar)? ff al. ^93 2 "En YL j ((Fo) ] (%o) 1) -20-Ar iss 


where LZ and [ Ay] have been defined in conjunction with equa- 
tion (A12). 


A similar process produces, for the right-boundary half- mesh, the 
formula | 


o B ï ho o “La Oh, (425) 
l" Gat mi A 


Now consider Ll case in which the applied moment is a step moment; 
that is, Mo(T) = 1. Equation (A24) reduces to 


0x 1 2053 E Dj 
A l a | (v "m 
1 3 " 


for left-boundary half-meshes. A similarity between the boundary 
formulas (A25) and (A26) and the interior formula (A12) is apparent; 
indeed, formula (A12) may be used everywhere if the sums of the values 
&t points 2 and 3 are simply replaced by twice the values at the inboard 
point of each boundary mesh. 
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For this problem, as in the preceding one, the magnitudes of the 
discontinuities & and SY at the wave front must be determined in 
advance for every point on the grid lines defining its position. From 
the boundary conditions and equations (18a) it is apparent that 
&»(0,0) = 8M(0,0) = 1 and $8V(0,0) = 0. Thus, for the line 7 - E, 
equations (20) become 


Sw cos AT 


(A27) 


Oo 
«il 
i 


-24 sin AT 


Furthermore, it is found that, when the wave front is reflected at 
either of the simply supported ends, the signs as well as the magnitudes 
of the discontinuities remain unchanged, in contrast to the behavior of 
a discontinuous wave front reflected at a free or fixed end as in the 
preceding case. Thus equation (A27) must be valid for the entire zigzag 
path defining the position of the wave front. 


With the discontinuities known, the step-by-step solution can be 
begun as before. Again, it is convenient to complete each triangular 
region before proceeding to the next. 


Computations for this problem were made on the Bell computer, with 
a 20-segment grid, for a beam with X= 5. The quantities @ and V 
were found in the first eight triangular regions; in addition, the moment 
at the center was computed. The time histories of shear at the end £= 


and of moment at the center £= 4 have been plotted in figures 5(a) 
and 5(b), respectively. 

Now consider the case in which a ramp-platform moment is applied to 
the beam; that is, Mo(T) is defined by 


Molt) = 7 (0S7É 
es (A28) 
Mo(t) = 1 (T >1) 
Equation (A24) reduces to 
04 1 em. CO) 1 
SARA "| LAMM | +2A (429) 
SE (à AT) “23 V. K Vy -2)FAr 
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for the left-boundary half-meshes for which T] 21 and reduces to 
eqvation (A26) for all the rest. Equations (A12) and (A25) are, of 
course, still valid for interior meshes and right-boundary half-meshes, 
respectively. 


Since 5w(0,0) = 8V(0,0) = O for this case, the solution begins 
with initial values o =V < O0 along the line r= E and there are no 
discontinuities to be added. Thus the solution may progress in any 
convenient manner without regard for the position of the wave front. 


The same quantities were computed for this problem that were obtained 
in the preceding case. As before, it was assumed that X = 5, and the cal- 


culations were performed on the Bell computer with a 20-segment grid. 
The resulting time histories of shear at & < 0 and moment at 6 = E 
are presented in figures 6(a) and 6(b), respectively. 
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- APPENDIX B 
` MODAL SOLUTIONS FOR UNIFORM BEAMS 


The solutions carried out in this and the following appendix make 
considerable use of Laplace transform techniques; all the Laplace 
transforms used, most of which were taken directly from references 15 
and 16, are given in table I for the sake of easy reference. 


Equations 


The equations of motion (2) may be written for & uniform beam 
vibrating in the absence of external distributed loading in nondimen- 
sionál form as 


e d We : EM aix 


Veg + (Ye - Y) = Vor = 0 


(81) 


where y - Y. The dimensionless bending moment M and vertical shear V 


D. 
are given by (see eqs. (1)) 


H= -V (Be) 
Y = R( Fe : y) (B3) 


If the beam is initially at rest (¥(£,0) = y (£,0) = y(E,0) = 
¥-(&,0) = O), equations (Bl) may be transformed to (where transform 6 in 
table I has been used) 


(BL) 
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oo 
-pT. 
in terms of the Laplace transforms Y(&,p) = J e Py(t,7)8T and 
O 


S . 
H(E,p) - L e 9Ty(5,T)à7. The solutions of equations (Bl) have the 
O E 


form 


Vip) cA T 


(85) 


V(t,p) = B(p)e™ 


Substituting these relations into equations (Bi) leads to a biquadratic 
equation in m 


"s i : EE (2) « 2) gE (B6) 


and gives the following relationship between B and A: 


WE "e 
ue 
B= ~ A (B7) 


Let the solutions of equation (B6) be written m = tim and im,. Then 


m 


mo 


and the general solution of equations (Bh) is 


Y(5,p) = C, cos më + Co sin mt + C4 cosh mos + Cy sinh m ( BB) 
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C4 M 
Y(£,p) = ^ m? + (=) (ce cos mE - Cy sin my E) + 


1 2 SAS 2 
tin mo^ - (=) p (Cy, cosh mË + C3 sinh mot) (B9) 


where C4, Co, C3; and C, are functions of p which must be chosen 
so that the boundary conditions are satisfied. 


Specific Problems 


Cantilever beam given a step velocity at the root.- If the 
C 


root £=0 of a cantilever beam is given a step velocity V = = Yr = 1 


at time 7 = 0, the boundary conditions may be written 


— 2 
y(0, T) x: €i 
w(0,T) = 0 
(B10) 
ve(1,7) = U 
y(1,7) = ye(1,T) 
These boundary conditions transform to (transform 8, table I) 

c2 1 
Y(0,p) = e E 
Vt (1,p) = 0 


g(1,p) = Y¿(1,p) 


The constants C4, Co, Cr, and Cj are determined by substituting 


equations (B8) and (B9) into equations (B11); the resulting expressions 
for Y(t,p) and ¥(&,p) may be written in the forms 
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“o Uy(£,p) 


bó sd d B12 
“1 p“D(p) Gs 


YE sp) 
and 


Co Ug(5,p) 


“1 p^D(p) mo 


Y(E,p) 


where 


œ | 
Uy(t,p) = Ç - E sin m; sinh m, + > cos mj cosh no) cos mié + 


LA 


mo 
í + m sin m4 sinh m4 + Q cos mj cosh m5jcosh Mo + 
l + my 3 
q sin m cosh m * m5 cos mj sinh mo sin m É - 


E sin mi cosh mo + Q cos HQ sinh ng) etm Mo É 


vg 
C 
Ug( E,p) = —— ll sin my cosh mo + Q cos mj sinh 02) (cos mé- 


cosh mot) - Q 5 Ll. 010 mj sinh mo + q COS M cosh mpjsin mjg + | 
l + = sin m4 sinh + Q cos m cosh sinh 
mj "ub E 1 mno mot 


ne = 


D(p) = 2 +- (a + à cos mj cosh m, + E - => sin my sinh m, 


e 2 
C 
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The inverse transforms F(&,7) and  w(&,T) are determined by 
substitution of Y(E,p) and y(£,p) into the complex inversion integral 


(ref. 15). The singularities of the resulting integrands eP'Y(E,p) 
and eP'y(t, p) must therefore be examined. Although the functions my 
and m, are, in themselves, multiple-valued functions of p with branch 


SUR EL 
points at p=0, EVB, and + — E, it follows from & considera- 
NC 
C 
tion of the fundamental theorem of the uniqueness of solutions of ordinary 
linear differential equations that Y(t,p) and v(f,p) must neverthe- 


less be single-valued. The integrands eP'Y(t,p) and e" y(E,p) thus 
have no singularities other than poles. The inverse transforms y(E,T) 
and wW(&,1t) will therefore be taken as the sum of the residues at the 


poles of eP'Y(t,p) and eP'g(t,p), respectively (see ref. 15). 


It can be shown that no singularities occur in the numerator 
functions Uy(&,p) and Ug(5,0); therefore, all the poles of Y(t,p) 
and (&,p) must be introduced by the zeros of the denominator peD(p). 
Consider first the equation 


s 
D(p) = 2 + (a + q)eos cosh + E - —|sin m sinh -0 (B14) 
This equation has an infinite number of roots p= ip,,: n-1,2,. .. 00 


all of which are, in general, simple poles of both Y(t,p) and qg(t,p). 
In addition, p =O can be shown to be a double pole of Y(£,p) but not 


a pole of Y(&,p). 


The sum of the residues of eP'Y(t,p) at the poles p= O and 
p= ipa, n21,2,3,.. . o, is 


= Da) (Es-Pn) - 
y(é,T) = T + » MP ie ePn™ , om =o) e ul 
nel 7 
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— similarly, the sum of the residues of eP'y(£,p) &t— p = 


tpn, 
= 1, 2, 3, .. . o, may be written 
. ) i Mp Pn) Po? , Up(5»-Pn) -pr 
WET) = Cy E Pr AM. dio e 
where 
: á 
D'(p) = = P] 
Rp =) 
Co 1 as ae 
mima (5 - alo m4 sinh Hs - 
— + hlcos m. cosh m; pa Ç z l + 
c 
m - m? : a M 
+ —— | >cos m sin mo - 
mo m? (m? N mo?) x 
2 
21“ 
SAL 
eo mo? = my " 
m ° (n, $ + mp?) mj? 
2 
Rp 2 
2 


— n] > sin m, cosh m 
2 2 l e 
no? (m T mo ) 
However, it is seen that Uy(E,p) = Uy(£,-p), Ug(t,p) = Ug(&,-»), and 
D'(p) - -D'(-p) so that these equations reduce to 
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E co Uy(&,Pn) 
y(E,T) = yr +2 ————— sinh p,T (B15) 
zl Pp D'(pp) 
and 
Co Uy(£,Pn) 
y(£,T) = 2 = — sinh p,T (B16) 


e 2 
1 n=1 Pn D' (Py) 


If, in equation (B14), p is replaced by ik, the result may be written 


dd r die icm (B1T) 


where 


and k= a, Q being the circular frequency of vibration. Equa- 

tion (B17) is the frequency equation for a uniform cantilever beam (see 
ref. 5) and its solutions k = ki, n= 1, 2, 3, « « . o, are the non- 
dimensional natural vibration frequencies of such a beam. Equations (B15) 
and (B16) may be written in terms of these natural frequencies simply by 
replacing p, by ik,. The results are 
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y(E,T) = al + 25 Xy, (6)sino er (B18) 


n=1 
and 
y(E,7T) = 22 y Xy, (£)s1n ky (B19) 
xb 
where 


Xy 20 mcm m = wa + = sin Pp sinh a e cos B4 cosh ag )cosh Ani + 
d - — = sin By sinh a, + Øn cos Bn cosh an)cos Bp - 
es cosh a, sin Bn + E sinh a4 cos Pn sinn an E + 


(= sinh an cos Bn + Øn cosh ap sin pr, sin Pnt 


se - (Bf 
"a Bn 

Xy, (5) = Enon a sinh a, cos By + Bn cosh An sin Bn (cos Bn& - 
Pn On : 

cosh Cta E) + Øn ah + En sin By sinh a, +t- 

1 Pn 

[m cos Bn cosh a, jetan Amé - Ç - on sin By sinh an + 


d. cos By cosh 2 d 
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2 
CL 
e 1 E + Bn 
Dn = —a (us - "d Bu sinh a, + 2 ae - hicos By cosh a, 
S , 
BnF 
2 
+ By®) 
NI S M aU. as 
2 2 
il j Rig E 
e(€1Y* 2/c1\? 
ihe zz) Bn? - Mn MAC A em B 
= = S apn COS B. 
an? (as? + By) m Bn (an? + B 2) 


ome (Rx 


The coefficients Xy. (E) and Xy (E) are seen to be the natural mode 


shapes associated with the frequency Kas As a check on the validity of 
the Laplace transform procedure these quantities may be shown to satisfy 
the differential equations (Bl) and the boundary conditions (B10). 


Substituting Y(€,7) and ¥(&,T) from equations (B18) and (B19) 
into equations (B2) and (B3) gives, for the moment and shear, 


M(£,T) = - I 3 Xy, (£)5in k,T (B20) 
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where 
s^ - (8) x? an 
Xy, (E) = CUR ZA + Ba sin Bn sinh An + 


B 
E cosh a, cos Bn) cosh ant - (2 - = sin B4 sinh a, + 


Dn cosh a, cos Bn) cos Bnk - Ex sinh a, cos B, + 


d. cosh ay, sin d (six Bak + = sinh ans) 
n 


V(£,T) = -2R z 1 Xy, (§) sin k,T (B21) 
L- E 


t - 2 sin B, sinh a, + d. cos By cosh cn} sn BR ~ 


P 
E sinh a, cos Bn + Øn cosh a, sin 3 (cos BE + i cosh as) 


The first eight terms in equations (B20) and (B21) have been used 
to compute the quantities M(O,T) and V(O,T). These results have been 
plotted in the range O $ r $ 10 in figures h(a) and h(b). 


Simply supported beam with an applied step end moment.- If a step 


moment M < L is applied to the end t=0 ofa simply supported beam 
at time T = 0, the boundary conditions are 
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y(0,T) = 0 
y(1,7) = 0 

(B22) 
yg (0,T) = -1 


Ye (1,7) = 0 


These equations transform to (transform 7, table I) 


(823) 


Ye (1,p) = 0 


If conditions (B23) are used to evaluate the functions Cj, Cos 
C4, and C, in equations (B8) and (B9), the quantities Y(£,p) 


and (&,p) become 


YLE ,p) = (cos m - cot mj sin mf - 
p(m? + mo 
cosh m5$ + coth mo sinh mañ (Bek ) 
2 
el 2 


GCE p) = - —— | sin mË + cot m cos me) + 
A m comm met) 


(sinh mot - coth mo cosh mot) (B25) 


42 | NACA TN 2874. 


or they may be written in the forms 


Ny(&,p) 


Y(€,p) = Dy) (B26) 


where 
Ny(&,p) = sin mj sinh mo (cos ma - cosh mot) - cos m, sinh mo sin m,& + 


sin My cosh mo sinh mot 


E 2 2 
Dylp) = p (my + Mo ) sin mj sinh m, 


and 
E Re, E,p) 
y(5,p) = Dy) | (B27) 
where 


N E,p) = mi E = da m (cosh m4 cosh nL - sinh mo sinh mot) - 


2 
a 
Co P 
2 SE m A : 
+= + 
mo m; =) 5 S no(cos m cos m$ sin m sin mE) 


Dg(p) = mymoDy (p) 


Here again, the functions Y(&,p) and U(E,p) are single-valued, 
end the inverse transforms y(E,T) and Y(t,T) of- the quantities Y(£,p) 
and T(€£,p) are taken as the sums of the residues at the poles of 
eP'y(t,p) and eP™Y(E,p), respectively. Consider all the roots of the 
denominators Dy Lp) and Dp(p). For mj = tnx or to = Tink, 


NACA TN 2874 43 


n=1,2, 3, . . . œ (any of which are roots of both denominators), 
equation (B6) becomes 


2 
do n e Gru 
eI C1 
[n m 1, 9. 3» x x 49) (B28) 


The solutions of equation (B28) are p = lia, and Tibus 
n = 1l, 2, 3, .. œ, where 


and these points are seen to be the locations of simple poles of both 
2 
eín = 


Y(t,p) and Y(€,p). The roots p < T E of the equation 


1 - [£2 
eL 

m? + mo^ =O are zeros of the denominator, but these are also roots 
of the numerators Ny(5,p) and Ny(f,p) and are not poles of either 
Y(5,p) or Y(E,p). The root p = O, on the other hand, is a pole of 
both Y(£,p) and Y(E,p) for, although both numerators also go to 
zero at this point, the denominators vanish more rapidly. The equation 
m4 =O has three roots: p = 0, which has already been discussed,. and 


p = HVR. The latter points are not poles of Y(t,p) since they are 
roots of Ny(&,p) = 0, but they are simple poles of {(£,p). 


First consider the residues of eP'Y(E,p) and eP'y(&,p) at the 
poles ‘tia, and tibpn. The sums of the residues at these poles provide 


the results 


| S ] N ib -ib a 
Yal 5,7) = J Ny( £, ian) lant da Ny(£, an) -i2nT " YLE n) ibgT Ny(£ n) -ibT 


ee ———— ———— : 829) 
LP m) O BpCim) DG) 5 By (“TBR ^ | 
and 
N (€,1ap) Lagt Ng(5, -iag) -lapnT Ng (5, ibn) ibgT Ny(é,-1by) ~IbyT ! 
Tute 3i IA 


where the subscripts on y and Y indicate that they constitute only parts - so-called "dynamic" 
parts - of the complete solutions for jy(£,T) and w(t,T). In order to obtain the values of 

N N 
x and. zx. at p = Tian and tib,, it is necessary to use the fact that, when p equals one 
Dy Dg! 


of the roots” tia, and Tib,, one of four equations m, = nx, My = -0N, My = ink, or m, = -inm 


is satisfied. The question of which value of p corresponds to each equation need not be answered 


since the following relationships can be derived: 


co 2 
(=) nr sin na 
Sy TR Ny c1 


(me) -= m em o uuu 7 
Dy' My =n Dy' m, --ns Dy moa lnn | Dy mj--inm- 2an b e * p? (an? + bn?) 


coy ene + aay nn § 
9m oo SCORES 
Py; M4 =DA i m4--ni 2) m-in1 D 


(B32) 
mi Pn Po + Pay? + Pye 


th 


«Log NL VOVN 
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As an example, let p= ia,- Then one of the four possibilities 
m, - nz, mj = D: m5- int, and m5 = -inn must occur. But for any 


Ny Ny 
of these possibilities, the ratios Dye and - have the forms given 


Dy 
by equations (B31) and (B32). Thus, 


coN2 
Ny(£, ian) E nx sin nn 
Dy'(1eg) l a L - bs?) 
(B33) 
Ng( 5, ian) E E M D 1| cos nx§ 
Dy' (ian) y an2 - bpe 


Ihe values of the other necessary ratios can be found in a similar manner 
and equations (B29) and (B30) reduce to 


co 


2 
a(t) = E) ya sin GF Sa _ cos =) (531) 
=j - anf - by An Ba 
and 
cos nxé by anf 
ya (6,7) = 2 a t= - co Ant - E - i) cos byt} (B35) 
n= 


Next, consider the contributions of the pole at p = O. Equa- 
tions (B24) and (B25) may be expanded in a Laurent series about p= O0. 
The resulting equations for Y(&,p) and W(5,p) may be written 


3 2 
Y(6,p) = je - E + = + vl (836) 


HEP) iE. ideo) (837) 


46 NACA TN 2874 


where O(p) signifies terms of order p or higher. Thus the quantities 
eP"Y(E,p) and eP'y(&,p) have simple poles at p=0 with the residues 


B 3 2 
CLE - +4 (B38) 
and 
y TN RS Ed (B39) 
st 2 2 3 R 


respectively. It can be readily verified that equations (B38) and (B39) 
actually constitute the solution of the problem when the beam is con- 
Sidered to be loaded statically. NE 


Lastly, consider the points p = tiVR, which are simple poles of 
V(£,p). It can be shown that 


(B40) 


so that the sum of the residues of eF'y(£,p) at-these poles provides 
the additional dynamic contribution 


Va (6T) =- ¿cos RT | (Bhi) 


Summing the contributions of ail the poles gives 


= 3 e CANE nx sin nn feos a T eog but 
SS — SM) (mo) 
2 3 C1 a = bp? &n by 
2 : p2 
WE) -i--teleL-lcosdET 759 met Ce ene an - 
) nzl an” = bn non 


zs (B43) 
eee Coe ae B43 
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Y ` 

The responses Jy(E,T) and Y(&,T) are seen to have the same form as 
solutions obtained by Williams' method (ref. 2); that is, they are the 
sum of a static part, Yst(E,7) and wYo+(&,T), and a dynamic correc- 
tion, Yal(£,7) and Ya, ( ET). E Yap ( £,T), which is expanded in a series 
of the natural mode shapes. The 'y- and Y-components of the natural 
modes have the shapes sin nx5& and cos nxë, respectively, and there 
are seen to be two natural frequencies, a, and bn, associated with 
each integer n so that each integer corresponds to two modes. From 
equations (B2) and (B3), 


oo 
2 2 
— nx Bin nxẸ || Dn an 
M(£,)-1-t&-e2 ze stn et o ajos T - - ljcos by 
s an - bp? nee "n nero 


(BÀ) 


and 


oo 


V(E,T) = cos YRT- 1 + OR E 2a coe &QT - COS byT) (B45) 


2 
an - By 
=] 


Six terms (n = 3) have been used in the summations of equations (B44) 
and (B45) in computing the quantities M(1/2,T) and “V(O,T). These 
results have been plotted in figures 5(a) and 5(b) up to += 8. 


Simply supported beam with an applied ramp-platform end moment.- 
Let the response to unit step moment at  & =.0, obtained above, be 
designated M4(5,7T) and V1( 5,7). Then the response to an arbitrary 


applied moment M(O ;T) may be obtained from Duhamel's superposition 
integral as 


M(&,7) = M(0,0)M,(&,7) + [ay E, T-@)Mg (0,0) de 
(B46) 


V(E,T) = M(0,0)V4(&,7) 2 V4(5,T-8)Mg(0,0)a0 
O 
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where the subscript @ indicates differentiation with respect to that 
variable. Let the applied moment have the ramp-platform time history 


M(O,T) = T (0S TS€1). 
= . (B47) 
M(O,T) = 1 (T>1 
Then equations (B46) become 
M(E,T) = L M,(€,7-0) a0 (o $T $1) 
sa (B48) 
W(t yr) - | M,(E,7-0)a0 (7 > 1) 
O 
and 
S | 
V(&,T) - | v, (8 ,T-0)a9 (0 S+Z 1) 
mE (B49) 
V(£,T) -| V4(&,T-8)ae (T 2 1) 
O 


Substituting M 
and V 


00 
2 . : 
= nz sin nz |/b sin ant 
Matti eae). A - E, 
nex 


and Ya from equations (B44) and (B45) gives for M 


n=l an? - bp? an 
2 in bnT 
ESA (0 Sr $1) 


* ee P , 
MET) =1-6 42 2 nr sin = | n em BT sin ag(T 1) E 
n= 


an? = bye nex? an 
an” sin bpr - sin bn(T - 1) 
Ss 2 Beb rc 20 ee (T > 1) 
none by 


(B50) 
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and. 
Ñ CO 
ES cos nz /sin aT sin byt 
Wt) > L sin Rr - Tea) Soa st (esp oot sip ba) 
(0 € 7T S 1) 


V(E,7) A A IR Tt - sin f - 1]- 2+ 


yR (B51) 


oo 
ER cos nng |sin ayT - sin a,(T - 1) 
n-I 8n? - bs? “n 


sin bpr - sin b,(T - 1) 
= E J (T >1) 


Bn 


As in the preceding case, computations using terms up to n= 3 in 
equations (B50) and (B51) have been made for M(1/2,7) and V(O,T) in 
the range O Z T $ 8. The results are plotted in figures 6(a) and 6(b). 
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APPENDIX C 


EXACT CLOSED SOLUTIONS FOR UNIFORM BEAMS WITH Cy = Co 


Introduction 


In this appendix some exact closed solutions of Timoshenko's equa- 
tions are derived for infinitely long uniform beams for which cy = Ga: 
The results for beams of infinite length are then utilized to obtain 
results for beams of finite length. This is done either by simply 
restricting the infinite-beam solution to-& time interval in which it 
coincides with the corresponding finite-beam solution or by superposing 
infinite-beam solutions in such a way that the result; satisfies the 
boundary conditions for a finite beam. 


Equations ` c 5 : 


If the propagation velocities cy and co are taken to be equal 
and if 'q-0, equations (14) reduce to 


De T M, = 0 (Cla) 
Ve IA. = 0 (C1b) 
EU TO. 
Mg +07 - V=0 (Clc) 
Ye - =V7 -5.=0 (cia) 
YS 

where XA = EE This system may be further XOUuce s. Dy MIXING & linear 

r 

i 


combination of equations (Cla) and (Clb) and a linear combination of 
(Clc) and (cla) - to obtain 


Zs + 21K, = 0 (C2a) 


Zy + 221Kg - 2112 = 0 (C2b) 


in terms of the complex variables Z = V + 20 and K =M + 9Aiv. 


"B 
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21 


For a beam initially at rest, equations (C2) may be transformed to 
(transform 6, table I) | 


T + 221ipk = O 


(C3a) 

Prin, + (p - 211) =0 (C3b) 
in terms of the Laplace transforms EL E.n) - | e P'z(t,Tr)àT and 
O 
k(E,p) = | e"PTx(E,mar. Eliminating Kk produces the equation 

O 

Cee - plp - 21) = 0 
which has the solution 
-EYp(p-211) \p(p-2A4) 
t(g,p) = &(p)e ^ + B(p)e* (ch) 


where A and B are governed by the boundary conditions. If the beam 
extends to infinity in the positive &-direction, boundary conditions 


stipulating that the beam remain undisturbed at E = 
Thus, equation (C4) reduces to 


eo must be satisfied. 


E(5,p) - Alp] EROM) 


(C5) 
Substituting this expression into equation (C3a) produces 
A -£ [p( 5-211) 
vs. = AP flp - ena) e n (c6) 


The quantity A(p) 


may now be determined from the remaining boundary 
conditions in conjunction with equation (C5) or (CÓ), or both. 


Specific Problems 


Infinite cantilever given a step velocity at the root.- If the 


end E =O of an infinitely long beam is restrained from rotating and 
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is given a unit step velocity at time T = 0, the boundary conditions 
are V(O,T) -.1 and o(0,T) = 0. In terms of Z and K these condi- 
tions may be written 


1.P.[Z(0,7)] = 0 
and 
I.P.[k(0,7)] = 2 


where I.P. designates the imaginary part of the quantity in brackets. 
The boundary conditions may be transformed to (transform T, table I) 


I.P. [£(0,p)] = 0 (CTa) 
and 3 

I.P.[&(0,5)]] = = | (CTb) 
respectively. l) 


Because of the nature of these conditions, it will be convenient 
to proceed in the following manner to determine the.function A(p). From 
equation (C5), it may be seen that A(p) = €(0,p); thus condition (CTa) 
establishes at once that A(p) for this problem is a real quantity. 
At E=0 equation (CÓ) becomes l ms 


or; in view of condition (CTb), 


R.P. [k(0, p} + e - AL fs (p E 2Xi) | (c8) 


where R.P. k(0,»)] is the real part of k(0,p). Since A(p) and p 
are réal, the conjugate equation, which must also hold, is 


R.P. [s(o, p)] - =< = «ee + 241) (C9) 
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Eliminating R.P.|k(0,p)] between equations (C8) and (C9) results in 


A(p) = au P=? B ER (C10) 


If A(p) is substituted into equation (C5), the result may be 
written in the form l 


-£ \p(p~2A4) 
C(&,p) =- lux? + eri(p - v»? * P) 
pip - 


which has the inverse transform (transforms 1, 3, 14, and 15, table I) 


z(&,7) = mete ola - EE) . 


"T å J [SXT - 0) 
d e P9; (192 - £8) (E — 1 a8 
E 
(tT >68) (C11) 


where Jo and J, are Bessel functions of the first kind. From equa- 
tion (C11) the shear and angular velocity are obtained as 


v(E,T) = -AF ¿cos - Jo(a 6? - t? ) + 


Jı BALT - 6 
[ sin A6 Jofa fo? e e = — 7 ag 
i (v > E) (c12) 


and 


@(t,7) = -2 {etm XT Joa 2 - £8) _ 
y LS Ji [Bx (T = 8)] 
Den ad AA as 


(Y > E) (C13) 
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in terms of integrals that apparently defy evaluation by analytical 
methods. i | id di 


The time history of the shear V(O,T) at the root of the beam 
may be obtained as a special case of equation (C12); however, it will 


be noted that Lív(o,7)) = £(0,p) = A(p) (where L{v(0,7)} denotes the 


Laplace transform of V(0,7)) so that the inverse of equation (C10) 
gives V(O,T) directly. Thus the shear at the root is found explicitly 
in terms of tabulated functions as (transform 13, table I) 


V(0,T) = “ine [cos AT Jg(XT) + sin AT J (AT) | (T >0) (Cib) 


Since A(p) may be replaced by L{¥(0 ,T)) , the inverse transform 


of equation (C6) may be written for this problem in the form (trans- 
forms 1, 2, 3, and 17, table 1) 


: 
K(&,r) = - $95 eth’ V(0,v-£) -| v(0,7-6)e ^7 


23 (1.10? - t? )lae (T >£) 


or, after integration by parts, 


" 
KET] = eriet®™ y (ar? - £2) ci fruc + 


EM sin A(T - 0) Jy [A(t - 2) Ss? - Glas 


(T >£) (C15) 


where V(O,T) is defined by equation (C14). From equation (C15) the 
moment and linear velocity may be written as 


-" T 
M(t,T) =--24 sin XT ov - e) - 2) i zx V(O,T-8)cos Ag - . 


Ed sin Afr - 6) Ji [Cr - oJ} JoG458 lag (r >e) (C16) 
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v(&,T) = cos AT 39 [M - 62) - IN à V(0,T-0)sin 19 + 


08 8 sin Mx - 0) zs [An - oro (de? -a Le >) (m 


again in terms of integrals requiring evaluation by approximate methods. 


As in the case of the shear, it will be more convenient to deter- 
mine the moment M(O,r) at the root by means other than a reduction of 
of the general relation, equation (C16). Since, for this problem, 


R.P. [«(0,»)] - L{M(0,7 J} , equation (C8) may be written, after substitu- 
tion for A(p), in the form 


L{M(0,7)} = BUS RAS ») 


The inverse transform of this equation is (transforms 4 and 14, table I) 


_ * 34 (208) 
M(O,T) = 22 —¿— 00 (Tr >0) (c18) 
O 


But, from reference 17, 


J1 (2X7) 


—— = 229, (2,1) - SIS 


and 


a 
a 
l Jg(238)a8 = son + = [31 (27) Hy (247) - oom n] 
where Ho and H] are Struve functions. Thus, equation (C18) becomes 


M(O,T) = SUE S nH1(227)| Jo(2ar) - [1 S xXrEo (2X )| n) 


(r»0) (C19) 
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Ihe Struve functions Hg and Hj are tabulated in reference 17 


in the range O S 2AT S 15.9. For larger values of the argument, it is 
convenient to use the approximations (ref. 17) 


Ho(2Xr) w Yo(2Xr) += —— (247 > 15.9) 
2 1 

Hy (2AT) s Y, (2x) + SU * ue) (24T > 15.9) 
A 


where Yo and Y, are Bessel functions of the second kind. From 
reference 18, 


Yo(9Xr)J4(2A7) - Y4(2Xr)Jo(2Xr) = — 


so that 


1 


pah a erp ST 


(SXT > 15.9) (C20) 


Computed values of V(O,T) and M(0,T), obteined from equa- 
tions (C14), (C19), and (C20) for a beam with X = 5, have been plotted 
in figures h(a) and (b), respectively, in the range O ST $2. In 
this range, the root of a finite cantilever beam behaves as if the beam 
were infinite, since the effect of the free end is ae felt until the 
return of the wave front at T = 2. 


Infinite simply supported beam with an applied step end moment.- 
If the end § = 0 of an infinitely long beam is simply supported and 


is subjected to & unit step bending moment at time T = Hun the boundary 
conditions are M(0,7) = 1 and v(O,rT) =0,’or K(0,T) 2 1. This 
condition transforms to (transform 7, table I) 


Fa 


kK(0,p) = = | = i (c21) 


MS) 


Substituting equation (C6) into condition (C21) reveals that 


21 


p Vp(p - 2i) 


A(p) (c22) 
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so that equations (C5) and (C6) become 


E(E,p) = SM nobis (ce3) 
dia Vs» - 211) i 


and. 


k(E,p) slet p(p-2M) (C24) 


respectively. The inverse transform of equation (C23) is (transforms 3 
and 15, table I) 


Bese). moie aci pre - 2) (T >£) (C25) 


and the inverse of equation (C24) may be written (transforms 2, 3, 4, 
and 16, table I) as 


2 2 

; E: ANO - & ) 

M oap | et i d 39 (v > &) (C26) 
t Ve? - £2 


From these equations, the shear, angular velocity, moment, and linear 
velocity are obtained as 


K(&,T) = et 


V(&,T) = -2X sin Xr I(r e. E) (T > E) (C27) 
O(&,T) = cos Xr Jol avr - x (T > E) (c28) 


M(t,T) - (t > E) (C29) 


> E) (C30) 


-] 


T 146? - z2) 
Y(&,T) = sin AB - A f sin AB do ( 
E 


ge - & 


Thus, V and are determined everywhere in closed form in terms of 
tabulated functions, but the relations for M and "v contain integrals 
which must be evaluated by approximate methods. 


58 NACA TN 2874. 


The results of this section are utilized in the next where a finite 
simply supported beam is considered. l 


Finite simply supported besm with an applied step end moment.- If 
the end 5 < 0 of a finite simply supported beam is subjected to & 
unit step bending moment at time T = O, the boundary conditions may be 
written i 


M(O,T) = 1 
EK = (C31) 
v(0,T) - M(1,T) - v(1,T) - 0 
or, in terms of K, 
K(0,T) = 1 
(C32) 
K(1,T) =0 E 


Although a direct solution of equations (C2) is possible in this case, 
the response of the finite beam may be obtained by a somewhat simpler 
procedure in which the responses of infinite beams to the same disturb- 
ance are superposed. This procedure is described below. 


A series of semi-infinite simply supported uniform beams extending 
in opposite directions are shown in figure 7. The beams have been posi- 
tioned in the figure so that the origins of the space coordinate E lie 
on the same vertical line, and the segments of greatest interest, 

OS Ets l, have been outlined with solid lines while the rest-of each 
beam is defined by dashed lines. For each beam the origin of the 
coordinate € is seen to lie at a different position relative to the 
end where a unit positive moment is suddenly applied at time T = O. 
An infinite number of these beams is assumed to.exist. 


The response of the top beam has already been determined as 
(eqs. (C25) and (C26)) 


Z(&,T) = exe 35 (7 fr - Z) (T >£) 


and 


K(&,7) = G(&, 7) (T > E) 
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where 


Tu e? - e? 
GL R.O = el^5 _ xk i Q1 31 We? - 57 ao 


g? . ¿l 


The response of the second beam may be determined from these equations 
simply by accounting for the changed direction and shifted origin of 
the coordinate Et. This process gives 


z(E,T) = ANE - (2 - "d (r»2- E) 


and 


K(&,T) = -G(2-€ ,T) (tT > 2 - E) 


In order to obtain the response of the third beam from the response of 
the first, only the shifted origin need be accounted for. Thus, for 
this beam, 


Z(&,r) = erie! ao ap” s(a 0? | (T >2 + E) 


K(€,7) = G(2+E,7) LT >2 + E) 


Similarly, for the fourth bean, 


Z(t,r) = zero | afr eoi "ui Lr = l 


K(£,T) = -G(4-E,7) (T >4 


E) 


E) 


and so forth, for all the rest. In each case, the region where the 
response is not specified is & region of zero response. 


Let the responses of all these beams to their respective dis- 
turbances be superposed. and consider only the response of the segment 
OSES1 of the resulting composite beam. Since the wave fronts of 

as 


all the disturbances travel at the same velocity > l, the response 


of this segment may be written 
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Z(&5,T) 


Z(E,7) 


Z(&,T) 


Z(E,7) 


and. 


K($,T) 
K( E, T) 
K(E,T) 


K(£,T) 


ote a - E2 ) 


210 T - £2) 4. 


c 
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LO ST S E) 
(R $T $2 - E) 
(2-&§S$TS2 + E) (033) 


221elAT nar - 2), Jo NE - (2 - "d + 


x 48 - c + 93] 


O 
a( E, T) 
G(E,T) a G(2-E,T) 


G(E,T) - G(2-E,7) + G(2+E,7) 


(2+-ESTS 


€ 
(2 - 8 


(24+ € 


A 
p 
i 


5) 


o 
^ 


S TSE) 


+ 
IA 


2 - £) 
2 + E) (C34) 


3 
HA 


TS - £) 


Since G(0,T) = 1, equations (C34) are seen to satisfy the boundary 
conditions (C32) and the response obtained above must be that of a 
simply supported uniform beam to a unit dimensionless step moment applied 
at the support € =O. 
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The vertical shear V and bending moment M may be obtained from 
equations (C33) and (C34) as 


I 


v(&,T) O LO $ T E) 


V(&,T) = -2X sin AT zb T 2) (EST 2 - E) 
V(e,r) = -2X sin XT tao S G 
I - (2 - 55 (2-8 


V(&,T) = -2X sin AT so? - Gl, ANE - (2 - "d 4 
«ME - @ v5] (9 +R ST S4- E) 


lA 
A 


2 + E) (035) 


T 


and 

M(&, 7) = 0 (OS TS E) 

M(&,T) = E(£,7) (EST $92 - E) 

M(E,7) = B(&,7) - H(2-E,7) (2-&S$7rS2+&) p(c36) 
M(t,7) = E(€,7) - E(2-E,7) + R(2«$,7) (2+ESTESÍ- $) 
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where 
p2 2 
JILAYOS - 
H(£,T) = cos AE - XE D cos A9 nil? - e) dae 
E Ve? - g? 


The angular and linear velocities œ and v may, of course, also be 
obtained from equations (C33) and (C34). 


Computations have been made for the quantity V(O,T) from equa- 


tion (C35) with X = 5. The results have been plotted in figure SL a) 
in the range OSr SB. 


Finite simply supported beam with an applied ramp-platform end 
moment.- Let the applied moment at the end E 


= 0 of a simply sup- 
ported uniform beam have the time history 


M(O,T) 


| 
-] 

7~ 
O 
IA 
+4 
IA 
H 

ene 


(C37) 
M(0,T) = 


| 
- 


(T > 1) 


The response of the beam to this disturbance may be obtained from the 
response to a unit step by using Duhamel's superposition integral. 
Thus, if the response to a unit step, as given by equation (C33), is 
designated Zj(E,T), the response Z(&,T) to the applied moment (C37) 
may be written 


Z(E,T) = IN Z4(5,T-0)d0 (o S T > 1) 


(C38) 
L 
Z(E,7) - 1 Z1(&,T-0)a80 (T > 1 
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Substituting from equations (C33) and letting & = O results in the 


following expressions for the end response 2Z(0,7T): 


Z(O, T) 


Z(0,7) 


Z(O, T) 


Z(0,T) 


Z(0,7) 


Ñ 
as || ei MT-8) [A(T - 6)|ae 


1 
= 221 l et 7-8) go lr - d ae 


1 
= 2. i T-6) 37 Tur - las 
i , [ar - Mao + 


-2 
MA i el Mr ol (T - 6)? - i lao 


1 
= 2M f eiM T-8) lar i ES + 


1 
hii l MESE 28)? = a 


1 
= 2M ix(T-0); |A(T - o)l ao 
l e JEG ) T 


. | 
USE) f SA Ege 1 + 


-4 
hi l EAR (t - 8)? - dj" 


(0S T€ 1) 


emm 
H 
HA 
-] 


IA 
-< 


(2 


HA 
+ 


(3 


A 
L 


(C39) 
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Tt can be seen that (transforms 3, 4, and 9, table I) 


"E: 
— Pri 
L | ego ( 18) as MENU ETT ES ¡== -1 
O pyp(p - 211) A dp P 


or (transforms 5 and 11, table I) 


œ | 
i el Jo (X0) a0 = Tel^T [Jo (7) - 137 (X7)| 
O 


Then, Z(0,T) may be reduced in the region O ZS T $2 to 


Z(0,T) = 2Xi L et 930 ( 28) dO 
Z(0,7) = 2rire WT atM) - 494 (27)] (0 $7 $1) (cho) 
and 
Z(0,T) = 21 f e Mia) as 
-1 


t | ST l 
211 i ei 39 (49) de -f eiM (10) 09 


Z(0,T) = exire!" [go (ir) -131 (47) | = 


SMT - dd Er DG - 1) - 191 [Cr - aj} 


(15T <2) (CH) 
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From these equations, the shear at £ =O may be written in this range 
as 


V(0,T) 


2Xxr|cos AT J4(ÀT) - sin AT atO (0S TS) 


v(0, 7) 


8X7 [cos AT 31 (A7) - sin AT Jo (Ar) | - 


(ch2) 
eX(T - 1)4cos AT - 224 [Mr - 2)] - 


sin A(T - 1)Jo NE - 3] (1 $7 $2) 


In the range T >2, however, the integrals in equations (C39) apparently 
cannot be evaluated analytically. If the shear V](0,T) due to a unit 


step moment has been computed, it may be convenient to write, for the 
shear at the point £=0 resulting from the disturbance (C37), 


S 
V(0,7T) - | V3(0,T -8)a8 (0 $ T € 1) 
O 


: (c43) 
V(0,T) - | v,(0,T-0)a8 (T > 1) 


1 
and then evaluate the integral l Vz(0 ,T-6)d48 numerically to obtain 


V(O,T) in the range T > 2. 


The quantity V(O,T), for a beam with X = 5, has been computed in 
the range OS T$ 2 from equations (Ch2) and has been obtained in the 


range 2 $738 by a numerical integration of the exact curve of 
figure 5(a) in accordance with equation (C43). These results are 
presented in figure 6(a). 
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TABLE I.- LAPLACE TRANSFORMS 


Number F(t) 
(a) 


ET - 6)Fp(6)a8 


Reference 


f(p) = E e"PTF(r)jdr 


£i(plfo(p) 


F(T - E) e^ SPf(p) 


el Mp(7) 


IN r(8)a6 (T » 0) 
0 


TR(T) 


f(p - 1X) 


= f(p) 


-£'(p) 


œd - pcd | (120) |zir() - Pireo) - ge 00) 2... 


1 


E (T >0) 
IQ(A7) 
J1(Ar) 
eir [ag (in) - 13,(17)] 
Dolar) + ul 


cos AT Jo(XT) + sin AT JT) 


2h 
— J) (247) 


5Y(T) =O for all values of T not specified. NACA 
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Figure 1.- Positive distortions and positive internal forces and moments 
associated with & typical besm element. 


X 


Figure 2.- The characteristics of Timoshenko's equations for a point 
in the x,t plane. 
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Figure 3.- Characteristic grid for application of nume 
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Figure 4.- Response of a uniform cantilever beam subjected to a step 
velocity at the root. 
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(b) Time bistory of the nondimensional moment at the root. 


Figure Àh.- Concluded. 
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(a) Time history of the nondimensional vertical shear 
at the end ¢£ = 0. 


Figure 5.- Response of a uniform simply Supported beam subjected to a 
step moment st E = 0. 


NACA TN 2874 


A G » — ——— — Numerico! solution (2O segments) 
O ( MOr — — — Modal solufion (static + 6 modes) 
t o | 2 3 
T 
e 
| 
Moment, M, O E 
at =$ 
-| 
Period of fist 
natural mode 
-2 
O | 2 3 4 5 6 ¢ 8 


Time, T 


(b) Time history of the nondimensional moment at the center. 


Figure 5.- Concluded. 
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Figure 6.- Response of p.uniform simply supported beam subjected to a 
ramp-platform moment at ¿£ =O, 
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Figure f.- Semi-infinite uniform beams superposed to obtain the response 
of a finite simply supported beam. 
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